Optical Activities as Computing Resources 
for Space-time Symmetries 



Y. S. Kinfl 

Center for Fundamental Physics, University of Maryland, College Park, 

Maryland 20742, U.S.A. 



Abstract 

It is known that optical activities can perform rotations. It is 
shown that the rotation, if modulated by attenuations, can perform 
symmetry operations of Wigner's little group which dictates the in- 
ternal space-time symmetries of elementary particles. 



1 Introduction 

Polarization optics serves as analog computers for the geometry of ellipsom- 
etry pQ and the Poincare sphere [2J. They also can perform algebras of 
two-by-two and four-by-four matrices known as Jones and Mueller matri- 
ces respectively [TJ [3]. It was established recently that these matrices cor- 
respond to the two-by-two and four-by-four representations of the Lorentz 
group which serves as the mathematical framework for Einstein's special rel- 
ativity USE]. 

Optical activities are known to speak the language of rotations. In the 
real world, all optical rays go through attenuations. If the attenuation is 
axially symmetric, it does not raise additional mathematical problems. 

In this paper, let us assume that the optical ray propagates along the z 
direction, and that the polarization rotates on the xy plane. If the attenu- 
ation along the x direction is different from that along the y direction, the 
combined effect of this asymmetric attenuation and the rotation around the 
z axis can perform an interesting mathematical operation. 
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It is a simple matter to construct a rotation matrix for a given value of 
propagation distance Z. So is the matrix for the asymmetric attenuation. 
However, the problem becomes nontrivial when these two effects are com- 
bined at a microscopic scale with a small value of z, and this combined effect 
is repeated to make up the finite value of z. 

It is shown in this paper that the resulting mathematics not only allows 
to make analytical calculations of the optical activities with asymmetric at- 
tenuation effects, but also provides a computational instrument for Wigner's 
little group which dictates the internal space-time symmetries of elementary 
particles. 

In 1939, Wigner noted a particle can have internal variables in addition 
to its energy and momentum [7]. For instance, an election can have its 
spin degrees of freedom, in addition to its momentum and energy. Photons 
can have helicity and gauge degrees of freedom. Wigner formulated this 
symmetry problem by introducing a three-parameter subgroup of the Lorentz 
group which preserves the four-momentum of a given particle. This subgroup 
is called Wigner's little group in the literature. 

For a massive particle, the little group is a Lorentz-boosted rotation 
group. For a massless particle, it is like (locally isomorphic to) the two- 
dimensional Euclidean group. For a tachyon with a space-like four-momentum, 
the little group is a Lorentz-boosted boost matrix, where the two boosts are 
made along perpendicular directions [7J [S] . 

It is now possible to understand optical activities in terms of Wigner's 
little group. Conversely, the optical activity can serve as an analog computer 
for internal space-time symmetries of elementary particles. 

In Sec. [21 we formulate the problem in terms of two-by-two matrices ap- 
plicable to the Jones vector. It is not difficult to write matrices performing 
rotations and attenuations separately. In Sec. [31 we compute the transfor- 
mation matrix if those two operations are performed at a microscopic scale, 
and are accumulated to a macroscopic scale. In Sec. [U it is shown that the 
transformation matrices correspond to those of Wigner's little group which 
dictates internal space-time symmetries of elementary particles. 
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2 Formulation of the Problem 



Let us start with a light wave taking the form 

(E x \ _ / Aexp{i(kz-ut + <f>i)}\ ( . 

\E y ) \Bexp{i(kz-ut + (P 2 )}J ■ { > 

This ray can go through rotations around the z axis. It can also go through 
xy asymmetric phase shifts and attenuations. The mathematics of these 
aspects is known as the Jones matrix formalism. If there are decoherence 
effects between the x and y components, the mathematics can be extended 
to the four-by-four Mueller matrix formalism [2]. 

It has been recently established that the Jones and Muller formalisms 
constitute the two-by-two and four-by-four representations of the Lorentz 
group [U [S] . In the two-by-two case, the transformation matrix is generated 
by three Pauli matrices 

* 1= (l J)' a2= (z 7)' a3= {l -l)' (2) 



plus three squeeze generators 




There are therefore six generators. They form a closed set of commutation 
relations. In mathematical terms, they form the Lie algebra for the SL(2,c) 
group which is locally isomorphic to the six-parameter Lorentz group. 

The three matrices of Eq. §Z§ generate the rotation subgroup of the Lorentz 
group. For possible computer applications, we are interested in the subgroup 
which produces real transformation matrices. Imaginary numbers are not 
convenient for computer mathematics. 

Among the above six generators given in Eqs.(T2]) and fl3j), cr 2 ,r 3 , and T\ 
are pure imaginary and can generate real transformation matrices. They 
satisfy the following closed set of commutation relations. 

[0"2,t 3 ] = 2zVi, [t 3 ,ti] = -2icr 2 , [ti,<t 2 ] = 2ir 3 . (4) 

This group generated by these three matrices is called the Sp(2) group and is 
applicable to many optical instruments and optical processes, either directly 



3 



or indirectly through its isomorphism with the S(l, 1) group for squeezed 
states. 

Within this framework, we are dealing with rotation matrices of the form 

R(6) = ( COSa - sin "V (5) 
V sin a cos a J 

and squeeze matrices of the form 

e x pH3))- «0 

and their multiplications. All the matrices in this representation are real. 

Indeed, optical activities can be described by these real matrices. The 
polarization goes through the rotation 

17 } Uin( 7 z) cos( 7 ^) J' 1 } 

as the ray propagates along the z direction. This matrix is applicable to the 
x and y components of the polarization, and the rotation angle increases as 
z increases. 

The optical ray is expected to be attenuated due to absorption by the 
medium. The attenuation coefficient in one transverse direction could be 
different from the coefficient along the other direction. Thus, if the rate of 
attenuation along the x direction is different from that along y axis, this 
asymmetric attenuation can be described by 

/ exp (-nxz) \ Xz f exp (fiz) \ , , 

I exp(-fi 2 z)J \ expt-^V' {Q) 

with 

A*2 + A*i 



A = 



2 



A*= o • ( 9 ) 



The exponential factor exp (—Xz) is for the overall attenuation, and the ma- 
trix 

' exp (fiz) 

exp (—fJ>z) 



(10) 



4 



performs a squeeze transformation. This matrix expands the x component 
of the polarization, while contracting the y component. We shall call this 
the squeeze along the x direction. 

The squeeze does not have to be along the x direction. It can be in the 
direction which makes an angle 9 with the x axis. The squeeze matrix then 
becomes 



S(9, fiz) 



cos 9 — sin 9 \ / exp (fiz) \ / cos 9 sin 9 

sin9 cos9 J \ exp (— jj,z) ) \ — sin# cos9 



which can be compressed to one matrix 

cosh(/i2;) + sin(2#) sinh(fiz) sin(2#) sinh(/z2:) 

sin(26 l ) smh(fiz) cosh(/iz) — sin(2#) sinh(/xz) 

If 9 = 45°, this matrix becomes 



(11) 
(12) 



S(n/A uz) = ( cosh ^ sinh M A (13) 

We shall work with this form of squeeze matrix in the following discussion, 
and use the notation S(ixz) without angle for the above expression. Thus, if 
the squeeze is made along the x axis, the squeeze matrix is 

S(0, uz) = #(-tt/4, /iz)S{nz)R(ic/4, uz). (14) 

If this squeeze is followed the rotation of Eq.(|7|), the net effect is 

cos(7^) — sm(jz) \ / cosh(/iz) smh(fiz) \ 



V sin(72;) cos^z) / \ sinh(/iz) cosh(/iz) / ' 

where z is in a macroscopic scale, perhaps measured in centimeters. However, 
this is not an accurate description of the optical process. 

This happens in a microscopic scale of z/N, and becomes accumulated 
into the macroscopic scale of z after the N repetitions, where N is a very 
large number. We are thus led to the transformation matrix of the form 

M( 7 , fi, z) = [e- Xz/N S{fiz/N)R{ 1Z /N)\ N . (16) 

In the limit of large N, this quantity becomes 

_ Xz/ \( 1 nz/N\( 1 - lz /Ny N 
\iizlN J \<yz/N 1 J. 



e 



(17) 
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Since '-jz/N and fiz/N are very small, 



-Xz 



1 
1 



M(7,/x,z) = 
For large N, we can write this matrix as 



_( 7 _^)x^ 

(7 + n) o ; N 



-,7V 



M(7, /i, 2) = e~ Xz exp (Gz) , 



with 



G 



-(7-/i) 



(18) 



(19) 



(20) 





( 7 + /i) 

The remaining problem is to calculate the exponential form exp (Gz) by 
making a Taylor expansion. We thus need to compute G N . This is a trivial 
problem if G is diagonal or can be diagonalized by a similarity transformation 
of a diagonal matrix. The problem arises because this is not always the case. 



3 Computation of the Transformation Ma- 
trix 

We are interested in computing the exponential form of Eq. ffTH]) . If 7 in 
Eq. (l20l is greater than /i, the off-diagonal elements have opposite signs, and 
we can write G as 

- exp (277) \ 

\exp(— 2r/) 

with 



G = Hexp(-2„) ")• < 21 > 



k = -v/7 2 -/j 
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or conversely 



exp (2 n ) = Jl±±, (22) 
V 7 - V 



7 = k cosh(2?7), /i = k sinh(2^). (23) 



If n is greater than 7, the off-diagonal elements have the same sign. We 
can then write G as 



6 



with 



k = y^ 2 -i 2 , 



(27) 



exp(277) = (25) 

or 

7 = fc sinh(2^), // = cosh(2^). (26) 

If 7 = mu, the upper-right element of the G matrix has to vanish, and it 
becomes 

(T 

,2 7 0, 

As /i becomes larger from fi < 7 to \l > 7, the G matrix has to go through 
this triangular form. 

We are now ready to compute the exponential form 

exp(zG). (28) 

The problem is whether it is possible to obtain an analytical expression of 
the above quantity. The usual procedure is to write a Taylor expansion. For 
this purpose, we need to calculate G N . We can manage this calculation when 
N = 2. However, for an arbitrary large integer N, it is not a trivial problem. 
This is exactly the problem we would like to address in this section. 
If 7 is greater than //, we write G of Eq. (I2T|) as 

with rj given in Eq. (l22|) . This is a similarity transformation of 



-1 

1 



(30) 



with respect to a squeeze matrix 



The role of this squeeze matrix is quite different from that of Eq. f|T3|) . It 
does not depend on z. 



7 



Let us go back to the G matrix, We can write G N as 











N 











and 



exp 



k 







cos(kz) 
sin(fcz) 



— sin(fcz) 
cos(kz) 



Thus, the exponential form exp (Gz) in of Eq. (l28|) becomes 











cos(fcz) — sin(fcz) 
sin(kz) cos(kz) 











and the transformation matrix of Eq. (fT6l) takes the form 



M( 7 ,/i,z) 



-Xz 



cos(kz) 
e~ 2v sin(kz) 



.-/It, 



sin(fcz) 



cos(kz) 



(32) 
(33) 

(34) 
(35) 



with k and rj given in Eq. 

If fi is greater than 7, the off-diagonal elements of Eq. fl20l) have the same 
sign, but we can go through a similar calculation. The result is 



M( 7 ,/i,z) 



-Xz 



cosh(&;,2) 
e _2r? sinh(fcz) 



e 2r] sinh(A;2;) 
cosh(fcz) 



(36) 



with k and 77 given in Eq. 

If 7 and fi are equal, the G matrix becomes 



G 



V2 7 



with the property 

/ 
V2 7 

and the transformation matrix becomes 



G 2 



(272 



(37) 



(38) 



(39) 



Let us go back to the case with 7 > fi. We can then gradually increase 
the parameter fi to a value greater than 7, which means from Eq. (1221) to 
Eq. fT25l) . This involves a singularity in the expression of exp (2r)) in these 
equations. This is a complicated mathematical issue [9J, but we can avoid 
the problem using the variables fi and 7. 
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4 Space-time Symmetries spoken by Optical 
Activities 



As was mentioned in Sec.flT]), the Lorentz group provides the basic mathemat- 
ical framework for polarization optics. The Lorentz group was used earlier 
by Einstein to formulate his special theory of relativity. In 1905, Einstein 
considered only point particles. After the formulation of quantum mechanics 
in 1927, it was found that particles can have internal space-time structures. 

If a given particle has internal space-time symmetries, such as electron 
spin and quark distribution inside a hadron, we have to rely on Wigner's 
little groups [7J. If the particle is massive, there is a Lorentz frame in which 
the particle is at rest. In this frame, the four-momentum remains invariant 
under rotations. However, its spin can change its orientation. Wigner's 
little group in this case is like (locally isomorphic to) the three-dimensional 
rotation group. We call this the 0(3)-like little group for massive particles. 

We do not observe particles with space space-like momentum or moving 
faster than light, but they play important roles in physical theories. We need 
those space-like particles in Feynman diagrams. For a particle of this type, 
there is the Lorentz frame where the momentum does not have its time-like 
component. It has its space component along a given direction. This four- 
momentum is invariant under Lorentz boots along the two perpendicular 
directions. The subgroup in this case is the Lorentz group applicable to two 
space coordinates and one time variable. We call this 0(2, l)-like subgroup. 

If the particle is massless, like photons, there are no Lorentz frames in 
which it is rest or with zero time-like component. For this case, Wigner in 
1939 observed that there is a three-parameter subgroup of the Lorentz group 
which leaves the four-momentum invariant, and that its algebraic property 
is the same as that of the two-dimensional Euclidean group. We call this 
£'(2)-like little group for massless particles. 

Let us now translate the two-by-two matrices given in Sees. [2] and [3] into 
the language of four-by-four Lorentz transformation matrices applicable to 
the Minkowski space of (z, y, z, t). In this convention, the momentum-energy 
four- vector is (p x ,p y ,p z , E) . If the particle moves along the z direction, this 
four- vector becomes 



in the unit system where c = 1, where m is the particle mass. We can obtain 




(40) 
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this four-vector by boosting a particle at rest with the four-momentum 



(0,0, 0,m), 



(41) 



with the boost matrix 



B(rj) 



/l \ 

10 

cosh(2?y) sinh(2?7) 

\0 sinh(2r]) cosh(2r/) / 



where 



tanh(2r^) 



P 



(42) 



(43) 



y/p' 2 + m' 2 

It is known that this boost matrix corresponds to the squeeze matrix B(rj) 
ofEq.(^U0. 

Now the four-momentum of Eq. (|4ip is invariant under the rotation matrix 



R(jz) 



( cos(2 7 2) sin(27^) 0\ 

10 

-sm(2jz) cos(27^) 

1/ 



Thus, the matrix 



B{Q)R(yz)B(-ri) 



(44) 



(45) 



leaves the four-momentum of Eq. (l40~l) invariant. While this matrix performs 
a rotation around the y axis in the particle's rest frame, we can also rotate 
this four-momentum around the z axis without changing it. This is what 
Wigner's little group is about for the particle with mass m. 

Although the matrix of Eq. (j4"5l) does not change the momentum, it rotates 
the spin direction of the particle in its rest frame. This is why the little group 
is not a trivial mathematical device. 

It is known that the rotation matrix of Eq. (T4"4"|) corresponds to the rotation 
matrix of Eq.([7]) [5j. Thus the two-by-two rotation matrix of Eq.flT]), together 
with the squeeze matrix of Eq. (l3T|) . generates the little group for particles 
with non-zero mass. 

If the particle has a space-like momentum, we can start with the four- 
momentum 

(0,0,p,£), (46) 
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where E is smaller than p, which it can be brought to the Lorentz frame 
where the four-vector becomes 



(0,0,p,0). 

The boost matrix takes the same form as Eq. (l42l . with 

. . E 
tanh(2^J = — . 

V 

The four-momentum of Eq. fj4"T|) is invariant under the boost 



(47) 



(48) 



S([iz) 



fcosh{2fiz) sinh(2/iz)\ 

10 

1 

V sinh(2/i^) cosh(2/iz) / 



(49) 



along the x direction. Here again the four- momentum of Eq. (14T|) is invariant 
under rotations around the z axis. 

The above four-by-four matrix corresponds to the two-by-two squeeze 
matrix of Eq.(fT3l) applicable to optical activities [5]. Thus, this squeeze 
matrix, together with the squeeze matrix of Eq. fl3T|) . generate the little group 
for particles with space-like momentum. 

Let us finally consider a massless particle with it four-momentum 



(0,0,p,p). 



(50) 



It is invariant under the rotation around the z axis. In addition, it is invariant 
under the transformation 



(51) 



This four-by-four matrix has a stormy history [8j [10] , but the bottom line is 
that it corresponds to the triangular matrix of Eq. fl3Tl) . and the variable 7 
performs gauge transformations. 

It is interesting to note that optical activities can act as computational 
devices for the internal space-time symmetries of elementary particles. 



/ 1 





-27 


2 7 \ 





1 








27 





1 -2 7 2 


2 7 2 


V27 





-2 7 2 


l + 2 7 2 / 
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Concluding Remarks 



Each human being has ten fingers. This is the origin of our decimal system. 
Vacuum tubes can do binary logic, and this is how electronic computer was 
developed. Quantum two-level systems can do more than the vacuum tube 
can. This is why we are interested in quantum computers these days. Indeed, 
computers are based on the mathematical language spoken by nature. 

Traditionally, polarization optics is known to produce the geometry of 
ellipse and that of the Poincare sphere. It also produces the algebra of two- 
by-two and four-by-four matrices. In this paper, we started with rotations 
combined with asymmetric attenuations in optical activities. It was shown 
in this paper that the optical activity speaks the mathematical language of 
Wigner's little group dictating internal space-time symmetries of elementary 
particles. 

We have used in this paper some mathematical methods not commonly 
seen in the conventional literature. In calculating the exponential form of a 
matrix, the usual procedure is to diagonalize the matrix by a unitary trans- 
formation. Then it is possible to write a Taylor expansion of the diagonal 
matrix. 

What should we do if the matrix cannot be diagonalized by a unitary 
transformation? Let us go back to the G matrix of Eq. fl2"Ul) . If 7 > /i, it was 
possible to bring the G the form 

where the similarity transformation matrix of Eq. (1311) is not unitary. It is a 
symmetric squeeze matrix. In addition, we used the property 



-1\ 2 (\ 0\ /0 -1\ 3 /0 -1 



1 J \o lj ' U J Vi 

to deal with the Taylor expansion. For fi > 7, we used 

1\ 2 fl 0\ /0 1\ 3 /0 1 



1 oj \o iy vi oj vi 

If fi = 7, the G matrix becomes triangular, and 



(53) 



(54) 



(T 2 
2 7 



0. (55) 
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The Taylor expansion truncates. 

Using these properties of two-by-two matrices, we were able to deal with 
the problem even though not all of them can be diagonalized. The triangular 
matrix of Eq. fl55|) is triangular and cannot be diagonalized. The matrix of 
Eq. fl54l can be diagonalized with the diagonal elements of 1 and —1. The 
two- by-two matrix of Eq. f[52"j) can also be diagonalized, but the eigenvalues 
are the imaginary numbers i and —i. However, the imaginary numbers are 
not too convenient for computer mathematics. Thus, we had to resort to the 
method presented in this paper. 
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